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Abstract 

Classical ghost imaging is a correlation-imaging technique in which the image of the object is found through intensity 
correlations of light. We analyze three different quality parameters, namely the visibility the signal-to-noise ratio (SNR), 
and the contrast-to-noise ratio (CNR), to assess the performance of double- and triple-intensity correlation-imaging setups. 
The source is a random partially polarized beam of light obeying Gaussian statistics and the image quality is evaluated as a 
function of the degree of polarization (DoP). We show that the visibility improves when the DoP and the order of imaging 
increase, while the SNR behaves oppositely. The CNR is for the most part independent of the DoP and the imaging order. 
The results are important for the development of new imaging devices using partially polarized light. 



Introduction 



r^Ghost imaging, or correlation imaging, is a novel imaging 
q technique in which information about the object is found 
"^through photon coincidences or classical light correlations 
Q-|[T]|2j . Originally demonstrated with entangled photon pairs 
^providing a high- visibility image |(3][4|) ghost imaging has 
since been realized with classically correlated [5] and inco- 
. ^herent (quasi- thermal) light |(6j[7)- Recently it has been ex- 
pended to imaging of temporal objects [8- 10 , of pure phase 

Essential 



^bbjects JTTJj an d through aberrations [|12 



16|. 



Q features of quantum ghost imaging can be emulated with 
L - — classical light, apart from the visibility as there is a constant 
^1 background present in classical imaging 



17 19]. The 



advantages of classical light are higher brightness and read- 
^^ily available sources. Various techniques have been intro- 
*sO duced to improve the image quality in ghost imaging. These 
include use of higher-order correlations [20 21], 
.ground subtraction | [22[ , differential ghost imaging 
"^■0 and computational ghost imaging 

The signal in ghost imaging is noisy and several pa 



back- 
[23], 



O 

(N 



rameters have been used to characterize the image qual- 



•ity. Two definitions are typically employed for visibility 



. [ |6|20 1 . Besides visibility, the signal-to-noise ratio (SNR) and 
^ the contrast-to-noise ratio (CNR) characterize image quality 
Hand various quantitative measures have been proposed in 
^ ghost imaging for both of them |22]|23j|26j|27j . While the 
visibility characterizes the contrast of the image, the SNR 
and CNR attempt to take into account the levels of the image 
noise. Most of the works concerning image quality in ghost 
or correlation imaging have dealt with scalar radiation and 
only recently the vectorial nature of light has been consid- 
ered by studying the effect that partial polarization has on 
ghost imaging with classical light [28-30]. In particular, it 
has been shown that the visibility increases monotonically 
with the degree of polarization (DoP) [ ]30[ and with the or- 
der of intensity correlations [28|, i.e., the number of arms 
in the ghost-imaging setup. 



In this work we investigate second-order (double- 
intensity, two arms) and third-order (triple-intensity, three 
arms) classical ghost imaging with partially polarized light 
obeying Gaussian statistics. Specifically, we analyze the im- 
age quality in terms of the visibility, SNR, and CNR. Both 
a general intensity-correlation arrangement and a specific 
ghost-imaging setup with lenses and an object in one arm 
are considered. We show that although the visibility im- 
proves with increasing DoP and is also better for third-order 
ghost imaging, the situation is not so straightforward with 
the SNR and CNR. On the contrary, due to increased noise, 
the SNR decreases with the DoP and the order of imaging, 
while the CNR is largely not affected by either. Owing to the 
noisy source, the SNR and CNR are generally small, even be- 
low one. 

The paper is organized as follows. In Sec. [2] we discuss 
electromagnetic intensity correlations and introduce the 
double- and triple-intensity ghost-imaging arrangements 
that we consider. In Sec. |3]we analyse the visibility in these 
setups and in Sees. |4]and|5]we calculate the corresponding 
SNR and CNR, respectively. A summary and a brief discus- 
sion of our results and their consequences are presented in 
Sec. [6] Some mathematical details are relegated to Appen- 
dices|AHc] 



2 Electromagnetic correlation imag- 
ing 

We begin by considering intensity correlation functions 
(ICFs) of random electromagnetic intensities correspond- 
ing to fields obeying joint Gaussian statistics. We then de- 
scribe the second-order (double-intensity) and third-order 
(triple-intensity) correlation-imaging setups. In both cases 
we first introduce a general intensity-correlation arrange- 
ment and then particularize it to ghost imaging with one 
object arm. 
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2.1 Degree of polarization 

Consider a stationary, uniformly polarized random beam of 
light propagating in the direction of the z-axis with the elec- 
tric field lying in the xy-plane. The beam's polarization state 
is characterized by the polarization matrix [ |31[ 
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where = (E*(t)Ej(t)), i,j € {x,y}, is the correlation 
function between the electric- field components and (...) de- 
notes the time (or ensemble) average. Due to the Hermitic- 
ity and non-negative definiteness of J there always exists a 
unitary transformation which diagonalizes J and the eigen- 
values Ji and J 2 are non-negative. Choosing J\ > J 2 the 
degree of polarization (DoP) is 



P 



Jl - J 2 
Jl + J 2 



(2) 



This expression is mathematically identical to the ones ob- 
tained by dividing the beam into fully polarized and unpo- 
larized parts and defining P as the ratio between the en- 
ergy of the fully polarized component and that of the total 
beam J31[ . Thus, when the light is completely polarized 
P = 1, unpolarized light has P = 0, and all other states are 
partially polarized. 

2.2 Intensity correlations 

The equal-time intensity correlation function (ICF) between 
N optical intensities is given by j3"Tj[32"|l 



,rjv) = (h 



IN 



(3) 



where I a = E* (r Q ,f) • F, a (r a ,t) with a G {1, ... , N} is the 
instantaneous intensity at point r a . For stationary light the 
equal-time ICFs are independent of time but depend on po- 
sition. The normalized version of G^ N ' is 
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(JV) 



(h) ■■■(!. 



(4) 



N 



which will be of convenience later. Each intensity can be 
expressed as I a = I ax +I ay = \E a}X (r a , t)\ 2 + \E a . y (v ai t)\ 2 , 
where E a ^ x (r a ,t) and E a ^ y (r a ,t) are, respectively, the x- 
and y-components of the electric field. Thus the Nth-order 
intensity correlation function can be divided into a sum of 



the Nth-order ICFs of the different components as [ 28 1 



(If- In) — 2^ (hh ■ ■ ■ Ini* 

ii,...,i N &{x,y} 



(5) 



where the sum is taken over i a G {x, y} for all a € 
{1, . . . , N}. Due to the underlying Gaussian statistics, we 
may apply the moment theorem and obtain ||20[|31|| 



(hn ■ ■ ■ Im N ) - ^2 r Mi " ' r 



nn_ 



(6) 
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where the summation is performed over all the N\ pos- 
sible permutations of the underlined indices and T^- = 
(E;,^,^^^,;)), with i,j £ {x,y} and a,j3 £ 
{1, . . . , N}, is the correlation function between the electric- 
field components. 



We can express the TVth-order ICF in terms of the equal- 
time mutual coherence (or the mutual intensity) matrix 



pa/3 



pa/? ~poL(3 

xx xy 

pa/3 

yx 1 yy 



(7) 



where the elements are defined as above. In the second- 
order we find that 



hh) = trr^trr 22 + tr(r 12 r^) 



(8) 



where tr denotes the trace and tiT a0! = (I a ) is the ath av- 
erage intensity. The first term is the product of the mean 
intensities (h) (h) and the second one is the correlation of 
the intensity fluctuations (Ah Ah), where AI a = I a — (I a ) 
with a G {1,2}. It is the latter term that is responsible for 
image formation in ghost imaging. Similarly, the third-order 
ICF takes on the form 



(hhh) = tiT n tiT 22 trr 



33 



trr dci tr(r^r' il ) 



+ ti-r 22 tr(r 13 r 31 ) + trr n tr(r 23 r 32 ) 

+ tr(r 12 r 23 r 31 ) + tr(r 13 r 32 r 21 ). (9) 

The first term is again the product of the mean intensi- 
ties (h) (h) (I3), the next three terms have a mean inten- 
sity multiplied by the correlation of intensity fluctuations 
between the remaining two intensities, e.g., (I 3 ) (AhAh), 
and the sum of the last two terms is equal to the correlation 
of all the intensity fluctuations (AhAhAI 3 ). 

2.3 Double-intensity correlation imaging 

A general second-order correlation imaging geometry is 
depicted in Fig. [if a). The source is a classical, quasi- 
monochromatic, uniformly partially polarized electromag- 
netic beam of light represented by the mutual intensity ma- 
trix r (r' 1 ,r 2 ) = J 7 (r' 1 , r' 2 ), where J is the polarization 
matrix and 7o(r' 1; r' 2 ) is a normalized spatial coherence func- 
tion. The source beam is split into two arms and propa- 
gated along the paths described by the kernels K\ and K 2 . 
The arms are taken to be polarization independent and their 
properties do not vary with time. The intensities h an d h at 
the end of the paths are measured as a function of position 
and correlated producing (hh)- 

The intensity correlation is given by Eq. ([8]). The mutual 
intensity matrices appearing in this equation can be written 
as ]29j||oJ 



JnT 



a/3 



(10) 



where 



70 (ri , r' 2 )K* (r Q , r\)Kp (r , r^dr^dr^ (11) 



with a, f3 G {1, 2}. Here and henceforth the primed coordi- 
nates refer to the source plane and the coordinates r a refer 
to the plane of the detector in arm a. We define the normal- 
ized form of T a/3 as 
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a/3 pa/3 



paap/3/3 



(12) 



a/3 I 



< 



which in view of the Schwarz inequality satisfies < |7 
1. Inserting Eq. Q10] ) into Eq. @ and employing Eqs. (|4f) and 
Q12| ) gives the normalized intensity correlation function 



5 (2) (ri,r 2 ) = l 



P 2 + l 

2~ 



(13) 
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where we used Eq. QA.4| ) of Appendix |X| The quantity 
has a constant background and a term that contains infor- 
mation on the correlation of the intensities I ± and I 2 and 
depends on the degree of polarization P. 
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Figure 1: (a) Generic and (b) specific double-intensity 
correlation-imaging setups. In (a), the source beam [r (r / 1 , r^)] 
is split into two arms characterized by the kernels K\ and Ki. The 
beam intensities Ii and J 2 at the end of the arms are measured and 
correlated ({hh})- In (b), a laser and a rotating ground glass disk 
(GD) are used to create a spatially incoherent light beam which 
then is divided into two arms with a beam splitter (BS). The ref- 
erence arm has a lens with the focal length / and the test arm 
contains the object with the transmission function T(r). The prop- 
agation distances are denoted by z iy i £ {a, b, c, d}. The reference 
arm has a CCD camera to measure the intensity distribution and 
the test arm includes a bucket detector (BD) which measures the 
total intensity. An intensity correlation ((/1J2)) is performed to 
form an image of the object. 



Double-intensity ghost-imaging setup We now consider 
the specific ghost-imaging setup shown in Fig. [TJb). The 
source is spatially completely incoherent, characterized by 
the mutual intensity matrix T (r[ , r' 2 ) = Jo^( r 2 — r i )> where 
S is the Dirac delta function. In addition, the reference arm 
(Kx) contains a lens and a high-resolution detector, and the 
test arm (A" 2 ) includes the object and a bucket detector with 
no spatial resolution. 

For a thin and large-aperture lens, the reference arm is 
described by the kernel J33fl 




(ri-ri) 



/ 2 . 2 
z b v\ + z a r{ 



(14) 



where z a and z b are the propagation distances, / is the lens 
focal length, ( = z a + z b — z a z b / f, and k = 2n/X is the 
wave number, with A being the wavelength. The transmis- 
sion function T(r) characterizes the object, which is pre- 
ceded with the propagation distance z c and succeeded with 



Zd- The kernel of the test arm is 

tf 2 (r 2 ,r' 2 ) 
1 



k 



Vz a z b J T(r)CXP l 2 



(r-r 2 ) 2 (r 2 - r) 2 

z c z d 




The normalized correlation between the output intensities 
is given by Eq. Q13] ), where 7 12 is specified by [see Eqs. Qll| ) 
and ([12])] 

f a/3 = j K* a {v a y)Kp{v p ,Y')dv' (16) 



with a,j3 € {1, 2}. If the imaging condition [34 35 1 
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holds, we find with the use of Eqs. fll4l)-([T7l) that 



(z a - z c y 

A Z b Z d 
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(17) 



(18) 



Thus an intensity correlation measurement can form an im- 
age (of the absolute value) of the object. The image is not 
due to the output intensities [trT 11 , trr 22 in Eq. (|8)] but 
rather results from the correlation of the intensity fluctu- 
ations at the outputs [tr(T 12 r 21 )]. Since < |7 12 | < 1, 
the normalized output intensity correlation is in the range 
1 < 9^ < (P 2 + 3) /2 suggesting that the image quality de- 
pends on the DoP. The range of holds both for the gen- 
eral geometry of Fig. [ija) and for the specific setup shown 
in Fig. [TJb) indicating that they have the same image qual- 
ity. 



2.4 Triple-intensity correlation imaging 

A general third- order correlation imaging geometry contain- 
ing three arms described by the kernels K\, K 2 , and K 3 is 
depicted in Fig. [2^ a). We make the same assumptions about 
the source (T ) and the arms as in the second-order case 
analyzed in the beginning of Sec. 2.3 The correlation of 



the intensities at the end of the arms, {IiI 2 I 3 ), is provided 
by Eq. ([9]). The mutual intensity matrices appearing in this 
equation satisfy Eqs. ([To) and ([TTJ with a,f3 e {1,2,3}. 
These equations together with Eqs. @, ®> ( [12] ), and QA.4| ) 
give the normalized triple-intensity ICF 



fl ( 3 )(r 1 ,r 2 ,r 3 ) = l 



P 2 + l 



(|7 12 | 2 + |7 13 | 2 + |7 23 | 2 ) 



3P 2 



(19) 



where . .) denotes the real part. We can express the 
complex correlation coefficients in terms of their magni- 
tudes and phases as 7"^ = |7 Q ^| e l *=' 3 . This implies that 

5J(7 12 7 23 7 31 ) = |7 12 ||7 23 ||7 13 | cos(A</> 123 ), where we have 
introduced the notation A</> 123 = <fri2 + 4> 2 3 — </>i 3 and used 
the fact that (7^) 



7 



P a . The function a' 3 ) has a con- 



stant background (equal to one) followed by three DoP- 
dependent terms which each contain information on the 



correlation of a different pair of fields [21 1. The last term 
has the strongest dependence on P and it carries informa- 
tion on all the three correlations. Because < |7 Q ^| < 1 for 



3 



all a,/3e {1, 2, 3}, the normalized third-order ICF is in the 
range 1 < < 3P 2 + 3 suggesting, as in the second-order 
case, that the image quality depends in an essectial way on 
the DoP of the source. 



ri = r 2 the triple-intensity ICF of Eq. ([19) then yields 



^(ri,n,r 3 ) 
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In other words, the information on the object is in the terms 

tr(r 13 r 31 ), tr(r 23 r 32 ), tr(r 12 r 23 r 31 ), and tr(r 13 r 32 r 21 ) 

correlator °f Eq. (|9), characterizing the correlations of the intensity 
. . fluctuations between the test arm and the reference arms. 

\J-lJ-2 3/ Because now 7 12 = 1 and < |7 13 | < 1, the normalized 

third-order intensity correlation of Eq. ( |21| ) is in the range 

(P 2 + 3) /2 < 5 < 3 > < 3P 2 + 3. 
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Figure 2: (a) Generic and (b) specific triple-intensity correlation 
imaging setups. Part (a) is similar to Fig. [TJ a) but has three arms 
instead of two. Likewise, part (b) is analogous to Fig. [ljb) with 
the exception that there is an additional reference arm which is 
identical to the first one. 



Triple-intensity setup with one object arm As the last 
arrangement, we consider the specific third-order ghost- 
imaging scheme shown in Fig. |2|b>). The source (r ), the 
reference arms \K\ and K 2 , as expressed in Eq. Q14D], 
and the test arm [K 3 , see Eq. Q15|)] are as in the specific 
double-intensity ghost-imaging setup of Se c. |2.3| The nor- 
malized ICF of the outputs is given by Eq. Q19|), where the 
functions 7°^ are obtained from Eqs. Q12] ) and ( [16] ) with 
a,/3e {1,2,3}. 

Since both reference arms are the same, the correlation 
function between them is 



il2 



exp 



ik 
~2 



1 - Zg/f 

c 



*(ri 



r 2 



(20) 



where the parameters are the same as those in Eq. ( |14| l. The 
Dirac delta function in this equation is due to the spatially 
fully incoherent source. Physically T 12 in Eq. Q20D repre- 
sents complete lack of correlation between the two fields 
when r! 7^ r 2 but shows complete correlation (with a uni- 
form amplitude of unity) when the points are equal. In 
Appendix |B] we use a spatially partially coherent source to 
explicitly demonstrate that the normalized correlation func- 
tion indeed is 7 12 = 1 when ri = r 2 . The remaining correla- 
tion quantities T 13 and T 23 are equal and given by Eq. fl!8| ) 
as long as the imaging condition of Eq. ( |17| ) holds. For 



3 Visibility 

Among the ghost-imaging quality parameters, the visibil- 
ity of the image has remained as an important quantity, 
although other parameters have been introduced. Visibil- 
ity (or image contrast) describes the relative difference be- 
tween the bright and dark areas of the image. Most studies 
on visibility in ghost imaging were initially performed in the 
double-intensity case Jl]|6]|19|l, but higher-order imaging 
has gained increased attention J20[|26|[36] | . Some prelimi- 
nary works on the influence of the degree of polarization on 
visibility have been carried out recently |28 |30[. 
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Figure 3: Qualitative sample realization of a normalized ghost- 
imaging signal. The object has two transparent regions with a nor- 
malized mean signal (? max corresponding to the bright areas of the 
image. The rest of the object is opaque leading to the signal (7 m i n 
and the dark areas in the image. The noise levels in the bright 
and dark areas are characterized by the root-mean-square values 
rms max and rms m i n of the fluctuations, respectively. 

Two main definitions are in use for the visibility in ghost 
imaging. The definition employed by Cao et al. [ j20[ for 
A^th-order imaging is 



V, 



where gmax and can be thought of as the average signal 
in the bright and dark areas of the ghost image, correspond- 
ing to the parts where the object is transparent or opaque, 
respectively. This idea is illustrated in Fig. [3] where a one- 
dimensional random realization of a signal related to an ob- 
ject with fully transparent and opaque regions is shown. 

The visibility employed by Gatti et al. |T] |19[| may 
be defined for second-order ghost imaging as Vq 2 ^ = 



(AO 



JN) 
ymax 



(jV) 
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(AJ 1 AJ 2 ) max /<7 1 J 2 ) 1 
ground, we obtain 



Since (Ji) (I 2 ) is a constant back- 
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There are two principal ways of extending this to higher 
orders, one of which assumes that (Ii ■ ■ ■ Pv) min = 
(ii) • • • (I N ) f$7\ and the other that 
(h ■ ■ ■ In-i) fe> (26). 



(If -In) van - 
Translated to the triple-intensity 

imaging schemes of Fig. [2] the first case corresponds to 
part (a); the minimum of the ghost-imaging signal is ob- 
tained when all intensity-fluctuation correlations disappear 
and only the first term in Eq. (J9) remains. The second 
case, (Iil 2 h) min = (hh) (h),^ closely related to Fig. |2^b), 
where the intensity correlation between the reference arms 
[the first two terms in Eq. ((9)] contributes to the back- 
ground, since these beams contain no information about the 
image. However, we want to take into account all the var- 
ious arrangements and thus we generalize the definition to 
the iVth order as 
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The visibilities are normalized so that < VX ,V, 



AN) 
ymax 

W rAN) 
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(24) 

< 1, 

but they scale differently between the end points. However, 
both definitions, Eqs. Q22| ) and ( |24| ), can be expected to lead 
to similar physical conclusions on the image quality. Next 
we make use of the results from Sec. [2] to calculate these 
quantities for double- and triple-intensity ghost imaging as 
a function of the degree of polarization. 

3.1 Double-intensity correlation imaging 

The correlation-imaging setups of Fig. [T] can both be de- 
scribed by the normalized second-order ICF of Eq. Q13| ) . Us- 
ing the facts that |7 12 | m in = and |7 12 | ma x = 1, the visibili- 
ties according to Eqs. ((22) and ((24]) then are |30j 
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v r 



(2) 



p2 
p2 
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(25) 
(26) 



Equations Q25| ) and ( |26| ) correspond to the maximum theo- 
retical visibilities and they are illustrated with the solid lines 
in Figs.[4fa) and[4fb), respectively. Both curves indicate that 
a rise in the degree of polarization affects the image visibil- 
ity positively. 

3.2 Triple-intensity correlation imaging 

The normalized output of a general triple-intensity 
correlation- imaging setup is given by Eq. Q191 ). Using the 
values | 7 Q/3 Uin = 0, |7 a/3 |max - 1 (for a,P€ {1,2,3}) and 
cos(A</>i 23 ) = 1, the visibilities from Eqs. ( [22] ) and ^§ take 
on their maximum values and are [ ]38 Q 



v c - 3p2 
"3P 2 



(27) 
(28) 



Equation ((27j) agrees with Eq. (17) in [J2Sjl and both 
Eqs. \27\ and ( |28D lead to the same physical conclusions 
that the visibility increases as a function of the DoP, as can 
be seen from the dashed lines in Fig. [4] The visibility in 
the general triple-intensity imaging scheme is considerably 
higher than that of double-intensity imaging, for each value 
of the degree of polarization. 
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Figure 4: Comparison of the maximum visibility according to 
the definition by (a) Cao et al. [Eq. ([22])] and (b) Gatti et al. 
[Eq. | |241 l]. The solid lines correspond to the double-intensity case 
[Eqs. ([25]) and ([26])], the dashed lines [Eqs. ([27]l and ([28])] and 
dash-dotted lines [Eqs. ( |29| l and < [30[ )] are for the general and spe- 
cific triple-intensity setups, respectively. 

Triple-intensity setup with one object arm Lastly, we ex- 
amine the visibility related to the specific arrangement in 
Fig. [2|b>) which has the normalized output of Eq. ( [21) . As 
shown in Sec. 2.4 in contrast to the general triple-intensity 
case, the normalized correlation between the first two arms 
in this setup is equal to unity (|7 12 | = 1) and it does not 
contribute to the image formation. For the remaining cor- 
relation parameters we have |7 13 | = | 7 23 |, | 7 13 | m m = 0, 
and | 7 13 |max = 1- Using these, the visibilities according to 
Eqs. ([22]) and ([24]) become 

5P 2 + 3 



V, 



(3) 



7P 2 

5P 2 



6P 2 



(29) 



(30) 



Equations ( |29| ) and Q30] ) are illustrated with the dash-dotted 
lines in Figs. |4]( a ) an d [4]M- The visibility of the spe- 
cific triple-intensity ghost-imaging setup behaves similarly 
to that of the two previous cases we considered above as 
the DoP is varied. It is smaller than the visibility in the gen- 
eral triple-intensity case but larger than in double-intensity 
ghost imaging, for all P. 

All the visibilities improve with the increase of the DoP. 
This is caused by the fact that while the backgrounds are in- 
dependent or weakly dependent on P, the correlation terms 
(which contain the image) show a stronger polarization de- 
pendence, as is evidenced for instance by Eqs. ( |13| ), ( |19| ), 
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and Q21] )]. For each P, the visibility in a given arrangement with this notation, 
is proportional to the relative abundance of the correlation- 
term contributions over the background. 

4 Signal-to-noise ratio 

The visibility is a measure of the (relative) difference be- 
tween the light (gmax) and dark (g^l) areas of the im- 
age. However, it contains no information about the noise, 
i.e., the strength of the intensity fluctuations within the im- 
age. We have assumed the source to exhibit electric field- 
vector fluctuations obeying Gaussian statistics. For such 
a source, ((AI) 2 )/(I) 2 = (P 2 + l)/2, showing that the 
higher P, the larger are the intensity fluctuations in rela- 
tion to the mean value |3T] 39 1. Hence the source noise 
level grows as P increases. The fluctuations are at their 
minimum ((A/) 2 ) = (I) 2 /2 when P = 0, indicating that 
the noise is comparable to the signal at the source. The 
source fluctuations will produce a certain amount of noise 
to the outputs, which we calculate for both double- and 
triple-intensity ghost imaging. Other sources of noise, such 
as those in the optical system or the detection, are not taken 
into account. 

In TVth-order ghost imaging, the quantity Ii ■ ■ ■ In is 
the fluctuating signal. Its variance is var [I\ ■ ■ ■ In) = 
{(h ■ ■ ■ In — {I\ ■ ■ ■ In)) 2 ) and the related noise is the 
square root of the variance (or the standard deviation), i.e., 
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(31) 



This is the root-mean-square (rms) of the deviation from the 
mean signal and it is always non-zero for Gaussian statistics. 
We define the signal-to-noise ratio (SNR) as 



snrW 



(/!•■■ I 



N 



I 



(32) 



N 



noise (I\ 

i.e., the average signal divided by the associated rms noise. 

The transparent parts of the object lead to bright areas 
in the ghost image with the normalized average signal g max 
and the average noise rms max , as is illustrated in Fig. [3] 
The maximum signal is obtained when the correlation pa- 
rameters in Eqs. Q13] ), Q19D , and Q21D take on their max- 
imum values. The SNR related to the bright areas is de- 
noted by SNR[ff^]. On the other hand, the opaque parts 
of the object lead to dark image areas, and these are char- 
acterized by g m [ n and rms m i n in Fig. [3] The signals for 
the dark areas are attained when the correlation parame- 
ters reach their minima, and the related SNR is denoted 
by SNR[ffj£2], In addition to the SNRs of the bright and 
dark areas, we examine the minimum and maximum SNRs 
within the ghost image, denoted by SNR^] and SNR^ , 
respectively. It can happen, for example, that the SNR maxi- 
mum is obtained when the signal achieves its minimum (i.e., 
SNRW «SNR[ffW]). 

For later use it is convenient, in view of Eq. Q31] ), to in- 
troduce the notation 



Figure 5: Dependence of the double- and triple-intensity SNRs 
on the degree of polarization. The solid lines correspond to the 
extremal SNRs of the second-order ghost-imaging setups shown 
in Fig. [I] The dashed lines depict the maximum and minimum 
SNRs of the general third-order correlation imaging arrangement 
[Fig. |2^a)], while the dash-dotted lines correspond to the ex- 
tremal SNR values in the specific triple-intensity ghost-imaging 
setup [Fig.|2fb)]. The maximum SNRs are plotted in purple and 
the minimum SNRs in green. 



4. 1 SNR in double-intensity imaging 

In Sec. H we calculated the second- and third-order ICFs. 
To find the SNR in the double-intensity case we also need 
the fourth-order ICF (l 2 I 2 ) or its normalized version g( 4 \ 
Since it is straightforward but tedious, the computation of 
the fourth- and higher-order intensity correlations is pre- 
sented in Appendix [c] For a uniformly polarized source 
and polarization-independent arms, an exact expression for 
the SNR as a function of the normalized correlation |7 12 | 
and the degree of polarization P is obtained by inserting 
Eqs. flT3] ) and QC4D into the definition given in Eq. fl34] ). 
From this expression we get the SNR in the dark (|7 12 | = 0) 
and light (|7 12 | = 1) areas of the image. The double- 
intensity SNR decreases monotonically as a function of \j 12 1, 
so it stays within the range 



P 2 



V5P 4 + 54P 2 + 21 



< SNR (2) < 



V-P 4 + 6P 2 + 5 



(35) 



n( 2N ) = 



{il...il)/(h) 2 ...(i N Y 



(33) 



where the upper limit corresponds to the dark areas 
and the lower limit to bright areas. Thus, for double- 
intensity ghost imaging we have SNR^ X = SNR[g^ n ] and 

which are shown with the solid 



(2) 



for the 27Vth-order ICF of N different pairs of intensities. 
The SNR can be written as 



snrW 



g(2N) 



[9 



(34) 



snr; 

lines in Figs. [5ja) and[5](b), respectively. The SNR is re- 
stricted in between these limits in the absence of other 
sources of noise. As opposed to the visibility, the SNR de- 
creases as a function of the DoR The mathematical reason 
for the SNR reduction with both the correlation parameter 
|7 12 | and the degree of polarization P is that, although the 



6 



signal increases with both |7 12 | and P, the level of noise 
grows proportionally more due to the stronger dependence 
of the higher-order intensity correlations on |7 12 | and P 
[cf. Eqs. fll3] l and flC.4] )]. Physically this is explained as 
follows: when P increases, the source noise grows as dis- 
cussed above, and when |7 12 | increases, the random output 
signals become more correlated and hence the fluctuations 
of their product grow. 



4.2 SNR in triple-intensity imaging 

Referring to the general triple-intensity ghost-imaging setup 
of Fig. a), we can obtain the associated SNR 



SNR (3) 



,(3) 



gW - U3)] 



(36) 



as a function of the magnitudes of the correlations between 
any two arms \j a0 \ (a, (3 e {1,2,3}, a ^ (3), the phase 
difference between the arms A^i 2 3, and the degree of po- 
larization P, using Eqs. Q19] ) and QC.5| ). 

First, we find the maximum and minimum signal-to- 
noise ratios SNR^ X and SNR^; n of the ghost image for 
each degree of polarization P. This is done computation- 
ally by maximizing and minimizing SNR (3) from Eq. 



at each P within the parameter ranges < \^ a ^\ < 1 and 
-1 < cos(A0i 23 ) < 1. Next the SNRs in the dark and bright 
areas of the image, SNR[gm2 x ] and SNR[g^jJj, are obtained 
with the assumptions |7 a/3 | m i n = and |7 Q ^| max = 1 and 
cos(A0i23) = 1, respectively. As in double-intensity cor- 
relation imaging, the dark regions of the image are found 



to have the maximum possible SNR, i.e., SNR[g^ n ] = 
SNR^ X , shown by the dashed line in Fig. [5ja) . However, 
in the general case the SNR related to the light regions 
(SNR^L]) is slightly higher than SNR^ n , illustrated by 
the dashed line in Fig. 5(b). In fact, SNR[gmax] of the gen- 
eral case equals SNR^; n of the specific case (discussed be- 
low), shown with the dash-dotted line in the same figure. 

The decrease of SNR^ when the correlation parame- 
ters and the DoP become larger is due to increased noise, 
as discussed in connection with the double-intensity case. 
Physically, SNR (3) is smaller than SNR (2) since the product 
of three partially correlated signals is noisier than the prod- 
uct of two such signals. 
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Figure 6: Dependence of the CNR on the degree of polarization. 
The solid line is for the double-intensity case [Eq. ( |39| l], while 
the dashed and dash-dotted lines correspond to the CNRs of the 
triple-intensity cases related to Fig. [5|a) [Eq. I |41| )] and Fig. [2jb) 
[Eq. I [42P ], respectively. 



5 Contrast-to-noise ratio 

The visibility measures the image contrast and the SNR as- 
sesses the noise in various parts of the ghost image. To study 
the contrast compared to the overall noise in the image we 
consider a third quality parameter, namely the contrast-to- 
noise ratio (CNR). It is defined as B22||27J 



cnrW = 



I 



N/ 



\ [ n0ise m, 



(h ■■■ I N ) + noise^ (Jj 



■In)] 
(37) 



Referring to Fig. [3j we see that the numerator in Eq. ( |3/p is 
related to the difference between the average signal in the 
bright and dark areas of the image and the denominator is 
the root-mean-squared average of the corresponding noises. 

5.1 CNR in double-intensity imaging 

The CNR of double-intensity ghost imaging becomes 



CNR (2) 



(2) _ (2) 
ffmax <?min 




=, (38) 



(2) 
,9m ax 



~(4) 
n ■ 
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where we have first substituted from Eq. Q31| ) and then nor- 
malized all the intensity correlation functions. Using the 
limits |7 12 | m in = and |7 12 | ma x = 1 for the normalized ICFs 
given in Eqs. fl!3| ) and flC.4] ), we obtain the expression 



CNR (2) = 



P 1 



V3P 4 + 30P 2 + 13' 



(39) 



Triple-intensity imaging with one object arm The SNR 

for the specific ghost-imaging setup of Fig. |2]^b) is obtained 
by inserting Eqs. Q21| ) and QC.6| ) into Eq. (|36|). Similarly 
to the second-order case, the SNR decreases monotonicaily 
as a function of |7 13 |, and thus its maximum and minimum 
values correspond to the SNR in the dark (|7 13 | = 0) and 
light (|7 13 | = 1) regions of the image, respectively. The 
SNR extrema are depicted in Fig. [5] by dash-dotted lines. 
We also found that the SNR in the bright areas of the spe- 
cific ghost-imaging setup is the same as that of the general 
triple-intensity correlation imaging case, while the SNR in 
the dark areas of the specific setup is lower than that of the 
general case. 



For fully polarized light the CNR achieves its maximum, 
CNR^ X = >/2/23 « 0.29. The CNR as a function of P 
is represented with the solid line in Fig. |6j It is almost con- 
stant with a slight rise for higher degrees of polarization. 



5.2 CNR in triple-intensity imaging 

The triple intensity CNR can be obtained analogously to the 
second-order case from Eq. Q37| ) and it is 



CNR (3) 



o (3) -o (3) 

i/max y min 
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For the general triple-intensity imaging case [Fi g. |2f a)], 
the normalized ICFs can be retrieved from Eqs. < \19y and 
QC.5] ). Then, inserting cos(A0 123 ) = 1 and |7 a ^| max = 1 

for ffl (3 ax and g£L, and | 7 Q/3 |mm = for g^ n and g^ n , with 
(a,/3 € {1,2,3}, a ^ the triple-intensity CNR with re- 
spect to the DoP is 



CNR (3) = 



12P 2 + 8 



V91P 6 + 1827P 4 + 3033P 2 + 577' 



(41) 



Equation ( |41D is shown with the dashed line in Fig. [6] For 
weakly polarized light the improved visibility in third-order 
correlation imaging leads to CNR values that are slightly 
higher as compared to the second-order case. For strongly 
polarized light the triple-intensity CNR falls below CNR (2) 
due to the increased noise. The maximum triple-intensity 
CNR is obtained with completely unpolarized light and it is 
CNR^ X = 8/V577 « 0.33. 

Triple-intensity setup with one object arm As for the 

general third-order setup, the CNR of the ghost imaging ap- 
paratus with two reference arms and one arm containing 
the object, depicted in Fig. [2|b>), is given by Eq. fl40| ). In- 
serting Eqs. ( |2ip and ( |C.6| ) related to the present setup, the 
CNR becomes 



CNR (3) = 



y/2(5P 2 + 3) 

V48P 6 + 947P 4 + 1602P 2 + 315 ' 



(42) 



where we have used the extrema I7 | mm = and 
I7 13 1 max = 1 for the correlation between the reference arms 
and the test arm. Equation Q42] ) is shown in Fig. [6] with 
the dash-dotted line. It behaves otherwise similarly to the 
general triple-intensity case [Eq. fl41] )] but is smaller for 
all P. This is because the contrast is lower (Fig. [4]) and 
the mean noise is slightly higher (Fig. [5]) for the specific 
setup than for the general case. The maximum CNR is 
CNR, (3 ^ X = y/2/35 ps 0.24, valid for completely unpolarized 
light. 



6 Summary and conclusions 

We have taken into account the electromagnetic nature 
of light and studied the influence the degree of polariza- 
tion (DoP) has on the performance of double- and triple- 
intensity correlation-imaging arrangements. The illumina- 
tion was a partially polarized, spatially incoherent beam of 
light obeying Gaussian field statistics. Both in the double- 
and triple-intensity setups a general intensity-correlation- 
imaging system and a specific ghost-imaging arrangement 
were considered. In all cases the image was assessed in 
terms of three different quality parameters, namely the visi- 
bility (contrast) and the signal-to-noise (SNR) and contrast- 
to-noise (CNR) ratios. 

We found that in all setups the image visibility improves 
when the degree of polarization increases and it is higher in 
third-order imaging than in the second-order case. In con- 
trast, the SNR was found to decrease with increasing degree 
of polarization as the source noise grows with the DoP. In 
addition, the SNR is lower for third-order imaging than for 
second-order imaging since in the triple-intensity case the 
image is created by a product of a larger number of noisy 
output signals. The CNR stays approximately the same as a 



function of the DoP, although it does go down slightly for 
the third-order imaging setups. If the CNR is considered as 
the most important quality parameter, as many authors do, 
then our results put into question the earlier suggestions of 
improved image quality with increased imaging order and 
DoP. Our results imply that increasing the imaging order 
might even be harmful and, on the other hand, changing 
the DoP does not have much of an effect and hence no po- 
larizers are necessarily needed in experiments. 

The SNR and CNR in ghost imaging with Gaussian light 
generally assume low values due to the noisy source. The 
impediments of the source noise are further enhanced by 
the correlations and the number of the output fields. The 
quality parameters studied in this work do not depend on 
the state of polarization of the polarized part of the illumi- 
nations, since the arms in the ghost-imaging setups were 
taken to be polarization independent. 

Acknowledgments 

This research was supported by the Academy of Finland 
(projects 128331 and 135030) and the Japan Society for 
the Promotion of Science (KAKENHI 23656054). 

A Traces of J 

In terms of the eigenvalues Ji and J 2 of J we have 

tr (J») = J? + J 2 " = J^l + x n ), (Al) 
(trJ )" = ( Ji + J2T = Jf(l + x) n , (A.2) 

where x = J2/J1 < 1. Using the degree of polarization 
[Eq. ((2)] we can write 



1 - P 

This implies that 

tr (Jg) (1 + P)" + (1 - Pf 
(trJ )" " 2" 

where we used Eqs. QA.1|)-(|aT3"]). 



(A3) 



(A.4) 



B Partially coherent source 

To calculate a value for the correlation coefficient 7 12 



in Sec. 2.4 we express the source coherence function as 



r (r' 1 ,r' 2 ) = Jo7(r' 2 — r^), where 7(r 2 — r[) is a normal- 
ized coherence function with the property 7(0) = 1. With 
the change of variables r = r 2 — r[ and R = (r^ + r 2 ) /2, 
Eq. ( |11| > becomes 

= J 7 (r)^;(r Q ,R-^)^(^,R+|)drdR. (B.l) 

Inserting the kernels from Eq. ( [14] ) corresponding to the ref- 
erence arms, the integral with respect to R is proportional 

to S [(rp — r a ) I (1 — Zb/f) — r] and we find that 



1 - Zb/f 

J (B.2) 




x exp 
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When ri = r 2 the normalized correlation takes on the value 

7 12 = f 12 /Vf n f 22 = 1. 



C Fourth- and higher-order intensity 
correlations 

To calculate the noise for the SNR and CNR in the double- 
intensity case we need to know the fourth-order ICF (I 2 I 2 }, 
and for triple-intensity imaging we need the sixth-order ICF 
(I 2 /|/|). These are obtained from Eqs. (J5j and ((6J) with the 
help of the matrix property 



E 

,n m £{x,y} 



A 1 A 2 ... A m — tr( A 1 

^ni» 2 H2"3 n n m n 1 ~ 11 V^ 1 



(C.l) 



where 



A* = 



A 1 . 



A 



A 



I 'I 

I 

yy 



are arbitrary 2x2 matrices for all z e {1, . 
the fourth-order ICF takes on the form 



(C.2) 

, m}. Therefore, 



</ 2 i 2 ) = [trr"] 2 [trr 22 ] 2 



tr(r n r n 



)[trr 



221 



tr^r 11 



)tr(r 22 r 22 



+ [trr 11 ] tr(r 22 r 2 
+ 2 [tr(r 12 r 21 )] 2 + 2tr(r 12 r 21 r 12 r 21 ) 
+ 4trr n trr 22 tr(r 12 r 21 ) + 4tr(r n r 12 r 21 )trr 22 
+ 4t r r n tr(r 12 r 22 r 21 ) + 4tr(r n r 12 r 22 r 21 ), 

(C.3) 

where T af} , with a,/3 e (1,2), are given by Eq. $7). When 
the source is uniformly polarized and the arms do not de- 
pend on polarization, the normalized version of Eq. QC.3| ) 
becomes 



duces to 



3P 6 + 39P 4 + 105P 2 



45 



(3P 6 - 
15P 6 



4 

111P 4 - 
- 231P 4 



225P 2 - 
- 285P 2 



45) | 7 13 
-45 



7 13 | 



(C.6) 



The generation of the sixth-order ICF and its subsequent 
simplification were handled with symbolic computational 
software. 
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